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Abstract 

A pragmatic approach for the construction of space-time codes over block fading channels is 
investigated. The approach consists in using common convolutional encoders and Viterbi decoders with 
suitable generators and rates, thus greatly simplifying the implementation of space-time codes. 

For the design of pragmatic space-time codes a methodology is proposed and applied, based on 
the extension of the concept of generalized transfer function for convolutional codes over block fading 
channels. Our search algorithm produces the convolutional encoder generators of pragmatic space-time 
codes for various number of states, number of antennas and fading rate. 

Finally it is shown that, for the investigated cases, the performance of pragmatic space-time codes 
is better than that of previously known space-time codes, confirming that they are a valuable choice in 
terms of both implementation complexity and performance. 

Index Terms 

Space-Time codes, block fading channels, performance evaluation, generalized transfer function. 

I. Introduction 

It is known since many years that the use of multiple receiving antennas, sufficiently spaced 
apart each other to obtain independent copies of the transmitted signal, is an efficient way to 
mitigate the effects of multipath propagation (see, e.g., [l]-[3]). However, only recently it has 
been realized that even the use of multiple transmitting antennas can give similar improvements 
[4]-[6]. With the introduction of space-time codes (STC) it has been shown how, with the 
use of proper trellis codes, multiple transmitting antennas can be exploited to improve system 
performance obtaining both diversity and coding gain, without sacrificing spectral efficiency 
[6]-[ll]. 

In particular, the design of STC over quasi-static flat fading (i.e., fading level constant over 
a frame and independent frame by frame) has been addressed in [8], where some handcrafted 
trellis codes for two transmitting antennas have been proposed. A number of extensions of this 
work have eventually appeared in the literature to design good codes for different scenarios, 
and STC with improved coding gain have been presented in [12]— [14]. In [15] it is pointed 
out that the diversity achievable by STC for binary phase shift keying (BPSK) and quaternary 
phase shift keying (QPSK) modulations can also be investigated by a binary design criteria, 
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instead of looking for the distances among complex transmitted sequences. This approach has 
been extended to multiple input multiple output (MIMO) block fading channels (BFC) in [16]. 

The determination of the STC with maximum diversity gain and largest coding gain remains 
a difficult task, especially for a large number of transmitting antennas and trellis code states. 
Moreover, the design of STC for fast fading channels is still an open problem. 

In this paper we present an approach to STC to simplify the encoder and decoder structures, 
that also allows a feasible method to search for good codes in BFC [17]— [19]. In fact, a criterion 
to achieve maximum diversity is given in [8], where, however, coding gain optimization is not 
addressed. Moreover, the STC in [8] require ad-hoc encoders and decoders. For these reasons, 
we present another possible approach to space-time coding, denominated pragmatic space-time 
codes (P-STC) [20]. Here, the "pragmatic" approach (the name following [21]) consists in the use 
of common convolutional encoders and Viterbi decoders over multiple transmitting and receiving 
antennas. We show that P-STC achieve maximum diversity and excellent performance, with no 
need of specific encoder or decoder different from those used for convolutional codes (CC); the 
Viterbi decoder requires only a simple modification in the metrics computation. 

We use the BFC model to investigate the design and the performance of STC. The BFC 
represents a simple and powerful model to include a variety of fading rates, from "fast" fading 
(i.e., ideal symbol interleaving) to quasi-static. 

Here, after the proposal of the P-STC structure, we first derive the pairwise error probability 
(PEP) of STC over block fading channels. Then, we propose a method based on suitable error 
trellis diagrams and generalized transfer function to evaluate a bound on the performance of 
STC over BFC, with a discussion on geometrical uniformity over the BFC. 

A new algorithm for searching good P-STC over BFC is then presented and applied to obtain 
the optimum (with respect to our performance bound) convolutional generators for various 
constraint lengths and fading rates. The numerical results, which compares our P-STC with 
the best known STC, confirm the validity of the approach. 

For simplicity we will focus on the BPSK and QPSK modulation formats, but the extension 
to other formats such as M-ary quadrature amplitude modulation (M-QAM) is straightforward. 

The paper is organized as follows: in section II the channel model and the general architecture 
of a system with STC are described; in section III the P-STC are presented; in section IV the 
PEP for STC over BFC is derived; in section V the frame error probability for STC over BFC 
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is analyzed; in section VI the search methodology for P-STC in BFC is illustrated; in section 
VII numerical results are provided, followed by the conclusions in section VIII. 

II. System architecture and channel model 

The general low-pass equivalent scheme for space time codes is depicted in Fig. 1, where n 
and m denote the number of transmitting and receiving antennas, respectively. 



We indicate 1 with C (t) 



Jt) J*) 

C x , . . . , On 



r 



a super-symbol, that is a vector of symbols si- 
multaneously transmitted at discrete time t on the n antennas, each having unitary norm and 
generated according to the modulation format by proper mapping. Thus, n symbols are sent in 
parallel on the n transmitting antennas. A codeword is a sequence c = (C^, . . . , C^) of N 
super-symbols generated by the encoder. 

This codeword c is first interleaved (we refer to intra-codeword interleaving) to obtain the 
sequence Cj = 1(c) = (C^, . . . , C* "^), where ai, . . . ,a N is a permutation of the integers 
1, . . . , N and X(-) is the interleaving function. 

The channel model includes additive white gaussian noise (AWGN) and multiplicative flat 
fading, with Rayleigh distributed amplitudes assumed constant over blocks of B consecutive 
transmitted space-time symbols and independent from block to block [17]— [19]. Perfect channel 
state information is assumed at the decoder. 

The transmitted super-symbol at time a t goes through the channel described by the (n x m) 
channel matrix = j/i^ j with i — 1, . . . , n\ s = 1, . . . , m, where h!f^ is the channel gain 
between transmitting antenna % and receiving antenna s at time a t . 

In the BFC model these channel matrices do not change for B consecutive transmissions, so 
that we actually have only L = N/B possible distinct channel matrix instances per codeword 2 . 
By denoting with Z = {Z 1; . . . , Z L } the set of L channel instances, we have 

H (CTt) = Z, for <r t = (l- 1)B + 1,...,IB , l = l,...,L. (1) 

When the fading block length, B, is equal to one, we have the ideally interleaved fading channel 
(i.e., independent fading levels from symbol to symbol), while for L = 1 we have the quasi-static 

'The superscripts H , T and * denote conjugation and transposition, transposition only, and conjugation only, respectively. 
2 For the sake of simplicity we assume that N and B are such that L is an integer. 
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fading channel (fading level constant over a codeword); by varying L we can describe channels 
with different correlation degrees [17]— [19]. 

At the receiving side the sequence of received signal vectors is r 7 = (R( CT1 \ . . . , R( CTJV )), and 
after de-interleaving we have r = J _1 (r 7 ) = (R^\ . . . , R^), where the received vector at 
time t is R^ 



rf • • • rm J with components 



n 

r? = V^EW+* 8 = l,...,m. (2) 

i=i 

In this equation is the signal-space representation of the signal received by antenna s at time 
t, the noise terms r$ are independent, identically distributed (i.i.d.) complex Gaussian random 
variables (r.v.s), with zero mean and variance N /2 per dimension, and the r.v.s hf], represent 
the de-interleaved complex Gaussian fading coefficients. Since we assume spatially uncorrelated 
channels, these are i.i.d. with zero mean and variance 1/2 per dimension, and, consequently, 
\hf].\ are Rayleigh distributed r.v.s with unitary power. The constellations are multiplied by a 
factor \/E~ s in order to have a transmitted energy per symbol equal to E s , which is also the 
average received symbol energy (per transmitting antenna) due to the normalization adopted on 
fading gains. 

The total energy transmitted per super-symbol is E ST = nE s and the energy transmitted per 
information bit is E b = E s /(hR) where h is the number of bits per modulation symbol and R 
is the code-rate. Thus, with ideal pulse shaping the spectral efficiency is nhR [bps /Hz]. 

For the discussion in the following sections it is worthwhile to recall that, over a Rayleigh 
fading channel, the system achieves a diversity V if the asymptotic error probability 
K( |M where K is a constant depending on the asymptotic coding gain [1], [22]. In other 
words, a system with diversity V is described by a curve of error probability with a slope 
approaching 10/D [dB/decade] for large signal-to-noise ratio (SNR). 

III. A Pragmatic approach to space-time codes 

In this section we present what we called pragmatic space-time codes, a low-complexity 
architecture for STC that allows an easy code design and optimization over fading channels 
[20]. The "pragmatic" approach consists in using common convolutional codes as space-time 
codes, with the architecture presented in Fig. 2. Here, k information bits are encoded by a 
convolutional encoder with rate k/(nh). The nh output bits are divided into n streams, one for 
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each transmitting antenna, of BPSK (h = 1) or QPSK (h = 2) symbols that are obtained from a 
natural (Gray) mapping of h bits. By natural mapping we mean that for BPSK an information 
bit b G {0, 1} is mapped into the antipodal symbol c = 2b — 1, giving c G {—1, +1}; for QPSK a 
pair of information bits a, b is mapped into a complex symbol c = (2a — 1)/ \[2 + j (2b — 1)/ a/2, 
giving c G {±l/-\/2 ± j/y/2}, with j = a/^T. Then, each stream of symbols is eventually 
interleaved 3 . 

We indicate the STC obtained with this scheme as (nh, k, fi) n-P-STC, where [i is the encoder 
constraint length and the associated trellis has N s = 2 fe(it ~ 1 ) states. For example, we report in 
Fig. 3 the four states (2,1,3) 2-P-STC encoder scheme for n = 2 transmitting antennas and BPSK 
modulation, obtained with a rate 1/2 convolutional encoder with generator polynomials (5, 7) 8 . 
We can describe P-STC by using the trellis of the encoder (the same as for the CC), labelling the 
generic branch from state Si to state Sj with the super-symbol C Sx ^ Sj — [ci, • • • ,c n ] T , where 
for BPSK q is the output of the I — th generator (in antipodal version). In Fig. 4 we report the 
trellis for the P-STC in Fig. 3. 

Similarly, in Fig. 5 we report the 4 states (4,2,2) 2-P-STC encoder scheme for n = 2 
transmitting antennas and QPSK modulation, obtained with a rate 2/4 convolutional encoder 
with generator polynomials (06, 13, 11, 16) 8 . 

It is clear now that with the pragmatic architecture the maximum likelihood (ML) decoder 
is the usual Viterbi decoder for the convolutional encoder adopted (same trellis), with a simple 
modification of the branch metrics. For example, in Fig. 6 we show the receiver architecture 
for the previous P-STC, that simply consists in the usual Viterbi decoder for the convolutional 
code adopted in transmission, with the only change that the metric on a generic trellis branch 
is J2T=i — V^l (hij s ci + h^gC-^j | 2 , being {q} the set of length n of the output symbols 
labelling the branch. In general, for n transmitting antennas, the branch metric for the Viterbi 
decoder is 

m n 

S = l 1 = 1 

Thus, we can resume the advantages of P-STC with respect to STC as in the following: 
• the encoder is a common convolutional encoder; 

3 In this paper we focus our attention on symbol interleaving: bit interleaving is addressed in [23]. 
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• the (Viterbi) decoder is the same as for a convolutional code, except for a change in the 
metric evaluation; 

• P-STC are easy to study and optimize, even over BFC. 
These aspects will be further investigated in the next sections. 

IV. The pairwise error probability for space-time codes over BFC 

In this section we address the performance analysis for the general class of STC over BFC. 
Given the transmitted codeword c, the PEP, that is the probability that the ML decoder chooses 
the codeword g ^ c, conditional to the set of fading levels Z, can be written as 



P {c - 9_\Z) = \ erfc^ ^-dP {c,g\Z) , (4) 

where erfc(x) = -7= e~* dt is the complementary Gaussian error function, and the conditional 
Euclidean squared distance at the channel output, d 2 {c,g\Z), is given by [8] 



n 
i=l 



cf- 9 f 



2 

(5) 



N rn 
t=l s=l 

To specialize this expression to the BFC we first rewrite the squared distance as follows 

TV m 

d 2 (c,g\Z) = (CW - GW) • (C« - G«) ff hf 

t=i 3=1 

N m 

^EE^^fe^r, (6) 

t=l s=l 

where = hi^h^, ...,/in,lj is the (1 x n) vector of the fading coefficients related to the 

receiving antenna s, and C^,G^ are the super-symbols at time t in the sequence c, and g, 
respectively. In (6) the (n x n) matrix 

A {t) (c,g) = (C« - G«) (CW - G^) H 

with elements 

4',!, = (<#'- 4°) W'-a?')* 

is Hermitian and non-negative definite 4 . 

4 This can be simply verified by noting that, since A can be written as A = yy H , for every (1 x n) vector x we have 
xAx H = xyy H x H = ||xy|| 2 > 0. 
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Due to the BFC assumption, for each frame and each receiving antenna the fading channel is 

s = 1, ... ,m, where is 



(!) „(2) 



,i L ) 



described by only L different vectors E j; 
the s-th row of Zj. By grouping these vectors, we can rewrite (6) as 



L m 



£ (c,g\Z)=J2Y, z s )F(l) ^9> 



(l)H 

S ! 



(7) 



1=1 s=l 



where 



FW(c >g )4 £ A<*>(c,g) = £ (CW-GW).(CW-GW) H Z = 1,...,L (8) 

teT(0 teT(0 

and T(l) = {t : = Z;} is the set of indexes t where the channel fading gain matrix is 

equal to Zj. This set depends on the interleaving strategy adopted. Note that in our scheme 
(Fig. 2) the interleaving is done "horizontally" for each transmitting antenna and that the set 
T(l) is independent on s, that means, in other words, that the interleaving rule is the same for 
all antennas. 

The matrix (c, g) is also Hermitian non-negative definite, being the sum of Hermitian non- 
negative definite matrices. It has, therefore, real non-negative eigenvalues. Moreover, it can be 
written as F^(c,g) = Jji 1 ) A^XJ^ H , where is a unitary matrix and is a real diagonal 
matrix, whose diagonal elements with i — 1, . . . , n are the eigenvalues of F^(c, g) counting 
multiplicity. Note that and its eigenvalues A® are a function of c — g. As a result, we can 
express the squared distance d 2 (c, g\Z) by utilizing the eigenvalues of F^(c,g) as follows: 



L m 



d 2 (c,g\Z) = ^^zfuC'AWU^zf" 

1=1 s=l 
L m 

= ^5]BfA«Bf ff 



where 



oil) oil) R ii) 
Pl^i > J 2,si ■■■■> yn,s 



1=1 s=l 
L m n 

= EEE A 

1=1 S=l 1=1 



(0 



(9) 



The difference between (9) and the similar expression reported in [8] is that, through (8), the 
eigenvalues in (9) are referred to the portions of the coded sequences with a given fading level. 
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Since represents a unitary transformation, has the same statistical description of z 



,(0 



Hence, in the case of Rayleigh distribution, has independent, complex Gaussian elements, 
with zero mean and variance 1/2 per dimension. Moreover, : 
independent V/ ^ j. Hence, the unconditional PEP becomes 



with zero mean and variance 1/2 per dimension. Moreover, for BFC, vectors and B^ are 



P {c -> g} = E < -erfc 



4iV 



m L n 



ft 



(10) 



s=l Z=l i=l 

where E {.} indicates expectation with respect to fading. By evaluating the asymptotic behavior 
for large SNR of (10) we obtain (see [24]) 



where 5 



P (c -> g) < K(mri) 



K{d) 



nn^ 

.1=1 i=i 



1 /2d - 1 



(0 



4iV n 



(ID 



2 M V d 

the integer rji = rji(c,g) is the number of non-zero eigenvalues of F^(c, g), and 77 (that we can 
call the pairwise transmit diversity) is the sum of the ranks of F®(c,g), i.e. 

L L 

V = V(c,g) = $>ank [F (l \c,g)} = ■ ( 12 ) 

/=i ;=i 
The PEP between c and (? shows a diversity mr? that is the product of transmit and receive 

diversity. 

Equation (11) can be seen as the generalization to BFC of the PEP for the quasi- static channel 
in [8]: for the BFC, to obtain the PEP we must compute the product and the number of non-zero 
eigenvalues of the set of suitably defined matrices F^(c, g) , I = 1, . . . , L, accounting, through 
(8), for the number of fading levels per codeword and for the interleaving rule. The analysis is 
valid for STC and will be applied also to P-STC. 

V. Error probability analysis for STC over BFC 

Given the transmitted codeword c the frame error probability, P w (c), can be bounded through 
the union bound as 

P w (c) <J2 V {^^9\ , (13) 



5 A looser bound can be obtained by observing that K(d) < 1/4. 
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that, using (11), gives for large SNR 



9¥=c 



E„ 



4iV n 



(14) 



.1=1 i=l 

where the dominant terms are those with minimum rj. It should be reminded that the parameters 
r/i and depend on codewords c and g. By retaining dominant terms only, the conditional 
asymptotic error probability bound becomes 



(c) « K(mr] min (c)) 



4iV 



-»?min (c)-m 



nn* 



(0 



(15) 



.1=1 i=i 

x} is the set of codeword 



E 

ge£(c i^min 

where r] min (c) = mm g r](c,g) and £(c,z) = {g^c:r](c,g) = 
sequences at minimum diversity. The asymptotic bound shows that the achievable diversity (also 
called diversity gain), i] min (c) -m, increases linearly with the number of receiving antenna. Note 
that, here, the transmit diversity order 77 min (c) has the same significant role of the code free 
distance, df, in AWGN channels. 

When dealing with codes for which the conditional error probability, P w (c), does not depend on 
the transmitted codeword c (see also the discussion in a following subsection), the unconditional 
error probability can be evaluated by arbitrarily selecting a reference codeword c . In the same 
way we may use P w (Qo) a s a bound for those codes for which we can prove that c is the worst 
case reference codeword. However, in general the error probability bound must be evaluated as 

^ = £ P {c} ^ (c) < £ £ P {c} P { c ^} , (16) 

where P {c} is the probability of transmitting the codeword c (i.e., for P-STC, equal to 2~ kN for 
equiprobable input bit sequence and 2~ fc ( Ar ~^ t+1 ) for a zero tailed code). By using (15), and by 
observing that the retained dominant terms are those with transmit diversity 77 min = minc^ m in(c), 
the asymptotic error probability bound can be written 



K(fj miQ m) 



4N 



-Vmin-m 



(17) 



E^fc} e nn* 

9G^(c,i?min) 11 = 1 1 = 1 

From (17) we observe that the asymptotic performance of STC over BFC depends on both the 



(I) 



achievable diversity, 77 min • m, and the performance factor 



^min(m) = £ P {c} F min (c, m) 4 P {c} Yl 

c c 9G£(c,r/ min ) 



nn* 

.1=1 i=i 



(18) 
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which is related to the coding gain in (17). 

Note also that 77 min and the weights Ylf =1 T\7=i \ f° r eac h Q. an ^ 9 do not depend on the 
number of receiving antennas. Therefore, when a code is found to reach the maximum diversity 
?7 m i n in a system with one receiving antenna, the same code reaches the maximum diversity 
?7min • rn when used with multiple receiving antennas. However, due to the presence of the 
exponent m in each term of the sum in (18), the best code (i.e., the code having the smallest 
performance factor) for a given number of antennas is not necessarily the best for a different 
number of receiving antennas. Thus, a search for optimum codes in terms of both diversity and 
performance factor must in principle be pursued for each m. 

To summarize, the derivation of the asymptotic behavior of a given STC with a given length 
requires computing the matrices F^(c, g) in (8) with their rank and product of non-zero eigen- 
values. In relation with [25], we also observe that: 

• By restricting in the bound the set of sequences g to those corresponding to paths in the 
trellis diagram of the code diverging only once from the path of codeword c, the union 
bound becomes tighter. 

• By restricting in the bound the set of sequences g to those corresponding to paths in the 
trellis diagram of the code diverging only once and only at time t from the path of codeword 
c, we obtain the first event error probability at time t. In the particular case of periodical 
interleaving over the BFC we can use the first event error probability at t = to obtain a 
simpler but looser union bound, in the following form: 6 

P w <Nj2^{c} J2 F i^^9}, (19) 

£ £e£o(c) 

where S (c) is the set of codewords g restricted to the first event error; this set must be 
used to evaluate the asymptotic performance (15) and the performance factor in (18). 

• From the error probability bound we can easily obtain an approximation by truncating 
the number of terms in the asymptotic expression (17) to the most significant terms, i.e, 
by keeping those terms with product of the non-zero eigenvalues smaller than a selected 
threshold 5p, or those terms corresponding to pairs (c, g) with Hamming distance smaller 
than a selected threshold 5h- 

6 This is also known as first event error probability analysis. 
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A. The New Concept of Space-Time Generalized Transfer Function for P-STC in BFC 

The evaluation of the error probability bound for P-STC can be carried out in an effective 
way, by extending the methodology in [24] for CC over BFC. This leads to the definition of 
the novel concept of space-time generalized transfer function (ST-GTF) for BFC. With respect 
to CC some modifications are required, as explained here, to account for the space-time fading 
channel. 

In order to define the ST-GTF let us first introduce the error sequences and discuss their role 
in the evaluation of error probability. P-STC are built using common binary convolutional codes, 
therefore they are group-trellis codes [26]. If we consider the input bit sequences and b g , of 
length kN, that generate the output codewords c and g, and define 7 e = b c @b g as the input 
error sequence for the transmitted codeword c and decoded codeword g, we can say that: 

- by encoding the input bit sequence e with the P-STC encoder we obtain a valid codeword; 

- given a transmitted sequence c (or the corresponding input b c ), the whole set of error sequences 
can be represented with the same trellis diagram used to describe the code. The all-zero path in 
this case describes the event of correct decoding. 

Having this in mind, we can rewrite the frame error probability bound as 

P w < Yl £ P P - C{b_ c © e)} , (20) 

where C(.) is the encoding function, P{& c } is the probability to encode the input bit sequence 
6 C , and P {C(b c ) — > C(6 C © e)} is the PEP related to input sequence 6 C and input error sequence 
e. As before, the bound is preserved by restricting the set of error sequences to those represented 
by paths in the trellis diverging only once from the all-zero path. The bound is also simplified 
by considering the error paths diverging from the all-zero path at t = 0. 

Thus, within this framework we can be proceed with the following steps to the definition and 
the exploitation of the ST-GTF of the code, for which an example is given in Appendix: 

a) construction of the error trellis diagram of the P-STC, starting from the trellis diagram 
of length N branches describing the P-STC. As observed, this trellis can be used both for the 
set of input sequences b^ and for the set of error sequences e (they both have the same trellis 
diagram of the convolutional code but with different meanings of input and output sequences). 

7 With © we denote the element-wise binary sum. 
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Let us denote with and the binary vectors representing the generic state at time t when 
the trellis is referred to the input sequence and to the error sequence, respectively. Recall that 
in our notation the number of states is N s = 2 k ^ l \ where /i is the constraint length of the 
convolutional code. Let us build the error trellis diagram according to [27, Chap. 12] by labeling 
the edges of the trellis referred to the error sequences e with N s x N s matrices ~E(se \ Se +1 ^), 
whose generic element i, j depends on the label Cs^Sj of the transition from state Si to state 
Sj, and is given by 



C, <?.._»,<?. — C . 



where C m „ (t+i) mC is the label of the transition from state si*- 1 © Si to si* +1 ^ © 

b) Construction of a modified error trellis diagram by labelling the generic transition —> 
Se +lS> of the error trellis with a new matrix label E'(si*\ Se +1 ^) whose generic element i,j is 
given by 8 

o-fc A (*) 

a] , (21) 

where Ai, . . . , Ajv are indeterminates and A^ is the n x n matrix given by 

A i% = (Cs^ Sj - C 8 w es ._ a (*+i) es .) (c Si ^ Si - C s («) es .^ a (*+i) es .) • 

This trellis diagram, named error ?re///s wzY/z error matrices, depends on the sequence of dummy 
variables A = (Ai, . . . , An) related to the multiple-input fading channel level as seen by each 
super symbol of the codeword. As before, due to finite interleaving, for each realization there are 
only a few different fading levels per frame. As an example, in case of quasi- static channels only 
one indeterminate must be used. In the opposite case of perfect symbol interleaving, although 
the number of indeterminates could be taken equal to the frame length, N, for the description 
of the average error probability over fading only one indeterminate may be used. 



c) Construction of the error trellis with error matrices for the BFC by using the same indeter- 
minate variable for super-symbols subjected to the same fading gain. This can be simply done 

8 In case of terminated codes by means of zero tailing the term 2~ fc has to be removed for t = N — fx + 1, N — 1. 
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with the position: 

A = J" 1 | A,..., A, . . .,p L ,.. : ,D L | , (22) 

B times B times 

which makes the trellis labels a function of the new set of dummy variables D = {Di, D L }, 
each related to one of the L fading levels. 



d) Evaluation, by using standard techniques, of the transfer function for the error trellis diagram 
as in [28], but with error matrices and by using the rules: 

Df' 1 ■ D^ 2 = £) Al+A2 (23) 
a-D^+b- L> Al = (a + b)D^ (24) 

a • A° = 1 (25) 

where A 1; A 2 are generic non-negative definite matrices, O is the zero matrix, a and b are scalars. 
In fact, we have to define for each node of the error trellis with error matrices, that is for each 
state = s at time t, a N s x 1 weighting vector polynomial Q S (D) which can be evaluated 
as the sum over all the transitions reaching s of the polynomials obtained by multiplying each 
transition label, which is a N s x N s matrix, by the weight of the node at time t — 1 from which 
the branch departs. Next, if we set to 1 the weight of the initial state of the trellis, denoted by 
O (the zero- state at the time 0), we can obtain what we call the space-time generalized transfer 
function (ST-GTF) as 

Tm(D) = Q Ojv (D) t U -1, (26) 

where Uq = [1 ... 0], O n is the final state of the trellis (the zero-state at time N) and the 
contribution of the correct sequence (the polynomial 1) is subtracted. 

In (26) the ST-GTF has the form of a polynomial in the indeterminates Di,...,D L with 
matrix exponents 

T M (D U D£) = W ( F(1) > ■ ■ ■ > F(L) ) • ^ W • • • D l (L) ^ 

(P(i),...,pW)#(0,... ) 0) 

where each pairwise error event is characterized by a set of L matrices (F^), . . . , F^) ^ 
(0, . . . , 0), and ^(F' 1 ', . . . , F^ L ^) enumerates (including the weight P {fe c }) the error sequences 
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producing F^, . . . , F^ L \ Among the terms in (27), the most important are those related to ma- 
trices F^, . . . , F^ L ) having minimum diversity, that is, the minimum value of Yla=i ran ^ [F®] . 
For these, it is important to evaluate the weight Ylf =1 nlli \^ given by the product of the all 
non-zero eigenvalues of F^, . . . , F^K 



e) Symbolic substitution of the powers in the ST-GTF with distances, by using the linear 
operator defined as 

L 



T 



a 



iia' 



F(0 



X (L) AL) 
a "1,1 a l,rn a L,l a L,7] L > 



(28) 



where a e M is an arbitrary number and A^, . . . , \$ are the rji non-zero eigenvalues of F®. 

With the same approach usually adopted for trellis codes, the ST-GTF in (27) can now be 
directly used to evaluate the error probability as 



P«, < 



^eJt[t m (d 



Li 



N 



(0|2 



(29) 



s=l 



This result is due to the well known bound erfc(x) < e~ x for x > 0. Tighter bounds and 
approximations can be obtained by using the results in [29], for example with the exponential 
bound erfc(x) < \(^ 2xl + \z~ x2 < e-~ x '\ or with the approximation erfc(x) ~ \ e ~ x ^ + \e~^ x2 . 
For large SNR the asymptotic union bound becomes, as in (17): 



K{mfj 



mm I ± mm 



m) 



E s 
4iV 



where 



F min (m)= ^(F (1) ,...,FW)- 

(F( 1 ),...,F< L ))G/ 



.1=1 1=1 



(I) 



(30) 



(31) 



Xf^ are the eigenvalues of F®, is the rank of F®, ^ min = min^^; and 

/={(F«,...,FW):X:^ = ^ 

is the set of error matrices giving rj min . 

We conclude the section with few remarks: 

1) The ST-GTF depends on both encoder and interleaver structures, which have been suitably 
considered to build the error state diagram specialized to BFC. 
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2) If we are interested in the evaluation of P w (c) conditioned to a selected reference codeword 
c (usually the one obtained with the all-zero input sequence) we need to define a ST-GTF 
referred to that sequence c , which can be easily obtained by using scalar (not matrix) labels in 
the error trellis diagram and the modified error trellis diagram. In this case, the generic transition 
s^P — > Se +lS> of the error trellis has to be labeled with Af- (t) where is the (n x n) matrix 
given by A^ = ( C (*) («+i) — C («),_ <t) (*+i)„ (*+i) ) ( C (*) (t+i) — C (t) m <*) (t+i)„ (t+i) ) 
and , . . . , s[ N ^ is the sequence of encoder states for sequence c (usually the all-zero state 
sequence). Moreover, at each state s^p = s the weighting polynomial is a scalar (not a vector), 
Q S (D), and the ST-GTF is simply obtained as T M (D) = Qo N (D) ~ 1. 

3) In a similar way, if the goal is to find the error probability according to (19) or to a tighter 
bound obtained by limiting the set of decoded sequences g to those corresponding to paths in 
the trellis diagram of code diverging only once from the path of codeword c, we can define 
a modified error trellis diagram by splitting the all-zero state at each time t, denoted by O t , 
into two states: O t , having only transitions departing to all the other states si* +1 ^ ^ O t+1 , and 
O t , having only the transition departing to O t +\ and all the transitions arriving from si <_1 \ By 
defining the time-t ST-GTF of this diagram as T Mt (D) = Qq jv (D) t U , when the initial settings 
are Q^^D) = (0, . . . , Of, Q 6t (D) = (1, . . . , 1) T and Q^D) = (0, . . . , Of for f + t, 
we can obtain: 

- the time — ST-GTF T M0 (D) whose use to evaluate (19) is straightforward, 

- the transfer function T' M (D) = J^t^o 1 ^Mt(D) which can be used in place of T M (D) to refine 
the error probability bound. 

An example of evaluation of the ST-GTF for P-STC over BFC is given in Appendix. 

B. Discussion on the geometrical uniformity for STC and P-STC 

Note that the error probability given in (14) is in general a function of the reference codeword 
c. The conditions under which there is no dependence on the transmitted codeword are related to 
the concept of geometrical uniformity, that has been introduced in [30] with respect to Euclidean 
distance. 

Geometrically uniform codes are codes with the same distance profile for all pairs of code- 
words. In AWGN channels, the geometrical uniformity guarantees that the performance is 
independent on the particular transmitted codeword. Thus, the frame error probability can be 
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evaluated by assuming the transmission of a particular codeword, that can be the 'all-zero' 
codeword generated when all the information (input) bits are 0. Clearly, this condition greatly 
simplifies the code design. 

However, the application of the concept of geometrical uniformity to STC requires a careful 
investigation, as highlighted in [31]. Indeed, it can be noticed that the PEP depends on the 
Euclidean distance of the coded signals after the multiple-input multiple-output channel an hence 
on the eigenvalues of matrices like those defined in (8), that can change with the reference 
codeword. For this reason, in general the design of STC should consider all possible transmitted 
codewords. 

For the P-STC introduced in section III we can easily see that: 



The P-STC before the channel (the set of codewords c) are geometrically uniform with re- 
spect to the Euclidean distance. In fact, for the P-STC with Gray mapping the Euclidean dis- 



tance between the symbols of two generic codewords, c, g, is o?e(c, g) = yJ2t I ? ~ dt^V = 
2^d H (c,g)for BPSK, and d E (c, g) = \/EtEi\c? ~ 9?\ 2 = V2^d H (^ , ~gj) + d H (c9 , ~gQ) 
for QPSK, where dn denotes Hamming distance and the superscripts /, Q refer to the real 
and imaginary parts, respectively. Since we are using convolutional codes, the Hamming 
distance spectrum is independent of the reference codeword, and therefore the same is true 
for the Euclidean distance spectrum. 

In a system with ideal symbol interleaving (BFC with L = N), the PEP in (10) depends 
only on the Hamming distance between the two codewords, but not on the specific reference 
codeword chosen. In fact, the PEP depends on the statistical distribution of the distance 
after the channel, defined in (5). In Rayleigh fading channels, each is a complex 
zero-mean Gaussian distributed r.v. with variance 1/2 per dimension; then the generic 
term hf\ {^f 1 — gf^j is still zero-mean complex Gaussian with variance 0.5 cf^ — . 
Note that the variance is thus proportional to the Hamming distance previously discussed 
between cf 1 and gf\ The resulting overall variable J2i h%s ~ 9?^ ^ s zero-mean 
complex Gaussian, with a variance that depends only on the Hamming distance between 
the codewords c and g. Since for ideal interleaving the r.v.s are i.i.d. also in t, we 
can conclude that the distribution of the r.v. defined in (5) depends only on the Hamming 
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distance between the codewords. 9 Thus, since for P-STC the Hamming distance spectrum 
is invariant with the reference codeword, the same applies to the error probability bound 
(13). 

• In the other cases and especially for quasi-static fading channels (L = 1), although the 
code is geometrically uniform before the channel, we could expect that the PEP depends in 
general on the reference codeword. In fact, it depends, through (9), on the eigenvalues of the 
matrix F^(c,g) defined in (8). However, in many cases we have numerically verified that 
the conditional error probability does not change significantly with the selected reference 
codeword. This happens in particular when: 

- The number of fading blocks L in the BFC is large enough with respect to the length of 
the error sequences; in this case the behavior of the ideally interleaved code is approached. 

- The memory of the code is small, and consequently the error sequences are short. In 
this case for many codes the distance in (9) has a distribution over the set of all possible 
codewords c which is mainly driven by the sum of the eigenvalues, i.e., the trace of F^)(c, g). 



This is again related to terms 



(t) (t) 2 
c\ -91 



and therefore to the Hamming distance between 
the codewords. Thus, the performance is mainly determined by the Hamming distance 
spectrum that, in P-STC, is invariant with the reference codeword. This will be verified 
numerically in section VII. 

Moreover, it is also worth noting that for P-STC with n = 2 antennas and BPSK modulation 
the error probability evaluated with the all-zero sequence as a reference codeword is always 
the worst-case error probability. 10 

In general we will not rely on the geometrically uniformity assumption (that holds before the 
channel but not after the channel), and so we analyze and design the P-STC by averaging over 
all possible transmitted sequences. We will also show, however, that fixing a particular reference 
codeword gives often similar results. 



9 It can be also shown that in (7) the matrix F ( ''(c, g) = A^(c,g) has only one non-zero eigenvalue given by X[ 1 ^ = 
J2i l c f — 9i^\ 2 directly related to the Hamming distance of supersymbols and G^. 

10 This can be proved (not included here for conciseness) by looking at the structure of matrix F^(c, g) and its eigenvalues. 
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VI. Search for the optimum P-STC on BFC 

In this section we address the issue of doing an efficient search of the optimum (in the sense 
defined later) generators for P-STC in BFC. Our search criterion is based on the asymptotic 
error probability in (17), so that the optimum code with fixed parameters (n, k, h, /x), among the 
set of non-catastrophic codes, is the code that 

- maximizes the achieved diversity, fj min ; 

- minimizes the performance factor F min (m); 

where the values of rj m - m and F min (m) can be extracted from the ST-GTF of the code. Therefore, 
an exhaustive search algorithm should evaluate the ST-GTF for each code of the set. 

Another search criterion for STC has been addressed in [12], [14] where a method based 
on the evaluation of the worst PEP was proposed. Although the worst PEP carries information 
about the achievable diversity, rj min , it is incomplete with respect to coding gain, thus producing a 
lower bound for the error probability. Even though our method based on the union bound is still 
approximate with respect to coding gain (giving an upper bound) it includes more information 
than the other method, leading often to the choice of codes with better performance. 

When applying our search criterion we must consider that, as shown in [32], the union bound 
for the average error probability is loose and in some cases (long codes and small diversity) is 
very far from the actual value. This problem can be partially overcome by truncating the sum to 
the most significant terms, but this technique leads to an approximation. However, this approach 
gives good results in reproducing the correct performance ranking of the codes among those 
achieving the same diversity f) min , as will be checked in the numerical results section. 

Of course, the achievable diversity is the most important design parameter. Since 77 min can not 
be larger than both n(c, g) < nL and the free distance df of the convolutional code used to build 
the P-STC, it appears that to capture the maximum diversity per receiving antenna offered by the 
channel, nL, the free distance of a good code for a given BFC should be at least nL or larger. On 
the other hand, there is a fundamental limit on the achievable diversity related to the Singleton 
bound for BFC [19]. In fact, if we define the reference block fading channel (RBFC) for the 
system as the ideal equivalent BFC with nL fading blocks that would describe the space-time 
fading channel if the n transmitters determine n independent channels, the achievable diversity, 
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which can not be larger than the diversity achievable on the reference BFC, is bounded by 



As an example, to achieve full diversity n with a P-STC on a quasi-static channel (L = 1) the 
value k/h can not be larger than I, thus the code rate of the convolutional code can not be larger 
than 1/n, or the value of h can not be smaller than k (see also [8]). 

Different methodologies can be used to compute the ST-GTF of the error trellis diagram: 
we can easily derive an error state diagram (by splitting all-zero error state) from the related 
trellis and in principle use classical techniques to evaluate T M o(D), but this approach is limited 
to long codewords with periodical interleaving since it could be computationally difficult to 
handle large matrices. The most efficient method to compute the ST-GTF is to proceed along 
the error trellis with an iterative algorithm which evaluates for each state s!P = s the weighting 
vector polynomials Q S (D) starting from t — 1 and ending in t = N with the initial conditions 
Qoq(D) = (1, . . . , 1) T and Q «»(D) = (0, . . . , 0) T . This method is also considered in [24]. 
Since 2 k N s branches connect the states of the trellis at each step, there are 2 k Ng products of 
polynomials (that can be reduced to 2 2k N 2 to account for zeros in matrix labels, and further 
reduced to 2 k N s if labels are scalars). 

With the view to utilize the ST-GTF to compare different codes in a systematic search for best 
codes, the previous algorithm still maintains a large complexity due to growth of the number of 
polynomial terms in the node weights when t increases, and only conventional simplification rules 
are available to reduce the evaluation complexity, which do not allow a significant improvement 
in the efficiency of the computation. 

This last issue is addressed in [24] for convolutional codes over BFC, where some simplifica- 
tion rules are given to largely reduce the computation complexity. Similar rules can be applied 
to derive the most significant terms of the ST-GTF, namely, those having small diversity order 
and product-degree, which allow the evaluation of ?7 min and F min (m). 

In order to formulate this method let us consider the error state diagram modified by splitting 
the all-zero states at each time t and, for each state s£^ = s, the weighting vector polynomials 
Q SjT .(D) which can be evaluated for the state s by using the initial settings Q s (o) ,^ o (D) = 
(0, ...,0) T , Qo t _ r (D) = (1,...,1) T and Q 6 (D) = (0, ...,0) T for t' ^ t - r. Let us also 
denote with W^(D) the set of N s vector polynomials Qs tjJ .(D), obtainable for each = s 
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at time t. The sets of vector polynomials Ur (D) can be computed along the error trellis, as 
for scalar weighting polynomials (see Appendix of [24]), with an iterative algorithm which 
starts from t — 1 and ends at t — N with the required initial setting for the all-zero state 
polynomial. We obtain at the end of the trellis T' M (D) = J2^=i Qd N r(-D) T U an d, if useful, 
Tmo(D) = Q 0jVjiV (Df Uo. 

The evaluation of the codeword error probability through the iterative computation for each r 
of the sequence Wr^(D), ...,Ur N \D) leads to the possibility of setting up much more efficient 
computation of a truncated asymptotic bound. To this aim we use the two following properties 
of non-negative definite Hermitian matrices [33]: 

PI) the rank of the sum of two non-negative definite Hermitian matrices is greater than, or 
equal to, the rank of each matrix; 

P2) the product of non-zero eigenvalues of the sum of two non-negative definite Hermitian 
matrices is greater than, or equal to, the product of non-zero eigenvalues of each matrix. 
Then, additional simplification rules are possible in order to eliminate polynomial terms which 
do not affect the final value of ?7 m i n and F min (m)). In fact, by means of PI and P2, respectively, 
at each step t it is possible: 

- to eliminate from each element of Q s r (D) the polynomial terms with rank of the exponent 
strictly greater than the minimum rank of the exponent of the polynomial terms in Q , r (D); 

- to eliminate from each element of vector Q s r (D) the polynomial terms with product of non- 
zero eigenvalues of the exponent much greater (a threshold should be fixed) than the minimum 
product of non-zero eigenvalues of the polynomial terms with minimum rank in Q , r (D)- 

VII. Numerical results 

This section presents the results of the previously proposed search algorithm used to design 
P-STC in BFC, and their performance (by simulation) in terms of frame error rate (FER) versus 
the SNR defined as E b /N per receiving antenna element. In addition, comparisons with the 
performance of previously known STCs are also given. All simulations are performed with 
random generation of information bits, thus without fixing a reference transmitted codeword, 
and with MIMO(n,m) we refer to a system with n transmit antennas and m receive antennas. 

First we investigate how the P-STC architecture exploits the diversity in BFC. To this aim, 
we evaluate the suitability of the pragmatic approach considering some P-STC obtained using 
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the known optimal convolutional code designed for the AWGN channel. As an example, for 
BPSK modulation with n transmitting antennas, a rate 1/n binary code is used, with a spectral 
efficiency of lbps /Hz. 

In Fig. 8 the FER for a BPSK modulated P-STC obtained with the de-facto standard 64 states 
convolutional encoder with octal generators (133, 171) 8 is shown, assuming a BFC. In particular, 
n = 2 transmitting antennas and m = 1,2, 4 receiving antennas are considered for various fading 
rates given by L = 1, 2, 5, 260 fading levels per codeword with N = 130. Note that even with 
this not-optimized choice of the generators, P-STC are able to reach the maximum achievable 
spatial diversity; in particular for L — 1 (i.e., quasi-static fading channel, meaning absence of 
time diversity), the diversity order is given by the product n ■ m. For L greater than 1, thus in 
the presence of available time-diversity, the achieved diversity order increases, depending also 
on the number of states. For MIMO(2,2) P-STC, typical values of interest for the FER (i.e., 
in the order of 10~ 2 ) can be reached with E b /N per receiving antenna element of about 6 dB 
(quasi-static channel), 3.8 dB (L = 2), 2.5 dB (L = 5) and 2 dB (fully-interleaved case), whereas 
with 4 receiving antennas the required SNR decreases to 0.2 dB, —1.1 dB, —1.6 dB and —2.1 
dB, respectively. 

The low complexity of the P-STC architecture makes also feasible the use of a larger number 
of transmitting antennas. As an example, the case of n = 4 is shown in Fig. 9 for quasi- 
static Rayleigh fading, iV = 130 and m = 1,2, A. Here, the convolutional encoder with optimal 
AWGN generators (135, 135, 147, 163) 8 is adopted [2]. Note that the case MIMO(4,2) achieves 
FER equal to 10~ 2 at 4.2 dB, that is greater than the 0.2 dB of the case MIMO(2,4) seen before; 
this is due to both the different power repartition on transmitting antennas and power combining 
at receiving antennas as well as the different code-rate. 

Similarly, by using the generators for CC over AWGN it is possible to design n-P-STC for 
QPSK (2 bps/Hz) by using the rate 2/2n convolutional codes. 

Let us now consider the search for optimum generators (in the sense defined in section VI). In 
Tab. I we report, for the quasi-static fading channel and QPSK, the characteristic parameters and 
performance of the best generators for the (4, 2, 2)2-P-STC, compared with the code proposed 
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in [8] 11 . Note that all codes achieve the available diversity n-m, but with different performance 
factors F min (m)/N. It is remarkable that the ratio between performance factors is almost the 
same as the ratio of the simulated FERs. Moreover, even at SNR of 15 dB, the asymptotic 
bound is sufficiently close to FER. Note also that the performance factor evaluated by fixing a 
reference codeword (F min (c , m) / N) provides a slightly different ranking of generators, giving 
as best code the generator (06, 13, 11, 16) 8 . This code is not the best according to F min (m) (that 
would give (05, 11, 06, 16) 8 ), but is very close to it terms of performance factor and the best in 
terms of FER. As will be clarified in the following, we checked that there are 11 codes out of 
2 16 behaving as the first and 47 behaving as the second, meaning that there is not a single best 
code but several codes that perform similarly. 

This fact suggests us to carefully investigate the performance differences among generators 
through exhaustive simulations. Thus, we performed an exhaustive simulation for all possible 
4 states n = 2, m = 2 P-STC in terms of FER for QPSK in quasi-static fading channel, with 
E b /N = 9 dB. In Fig. 10 we report the FER for all 4-states P-STC obtained through 2/4 
convolutional encoders (i.e., 2 16 generators that are ordered in abscissa). A remarkable outcome 
is that also for P-STC it is verified a phenomenon similar to what already discussed in [24] for 
convolutional codes in BFC: non-catastrophic codes can be divided in few classes, with almost 
the same performance for codes in the same class. Note that within the class of codes providing 
the best performance, there is the one obtained through our searching methodology, that gives 
a FER of about 0.01. Even for this simple case of 4 states generators, the exhaustive search by 
simulation required one entire week on a Pentium 4 personal computer, whereas with our code 
searching algorithm we saved about two order of magnitude in time. An exhaustive search for a 
larger number of states is impractical, while our search algorithm works still well, emphasizing 
the importance of algorithmic methods. 

Hence, it is important to note that the pragmatic structure is not only interesting from the 
implementation point of view, but it also provides interesting performance, that, in all cases 
we investigated, outperformed the previously known STCs. In order to make the comparison 
between P-STC and STC possible, in the following numerical results we assume iV = 130 [8], 

n It is possible to show that for these parameters the code given in [8] is amenable of a P-STC representation with generators 
(01,02,04, 10) 8 . 
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[12], [14]. As an example, the performance of our best MIMO(2,2) QPSK (4,2,2) 2-P-STC (with 
generators (06, 13, 11, 16) 8 ) is compared in Fig. 11 with other STC known in the literature for 
quasi-static channel [8], [12], [14]. These results show that out P-STC outperform previously 
known STC. 

Moreover, the proposed code search methodology also enables to find P-STC for various fading 
rates, that is, for various values of the parameter L. As an example, in the case of MIMO(2,2) 
QPSK (4,2,2) 2-P-STC, we obtain that the best generator is (05, 06, 13, 17) 8 when L > 2, as 
shown in Fig. 12 for L — 1,2, 5, and 260. At the author knowledge no other results for the BFC 
are present in the literature, so we compare our codes with the original STC in [8] even if the 
latter was designed for the quasi-static case. Note that with only 4 states codes are not able to 
exploit all available time diversity, but the proposed codes already achieve the available spatial 
diversity. 

Then, we investigate the impact of the number of states on the performance of MIMO(2,l) 
P-STC with BPSK in quasi-static fading channel. In Tab. II we report best codes obtained 
through the search algorithm for 2,4,8,16 states, for which we indicate the achieved diversity, 
the performance factor and the FER. We also report the performance factor for AWGN optimal 
generators with the same number of states. Note that all codes achieve the maximum diversity, 
and that increasing the number of states does not produce relevant performance improvements. 
Moreover, on the quasi-static channel the error probability bound tends to become looser, 
especially when the free distance of the convolutional code increases with respect to the achieved 
diversity. 

The behavior is different in BFC with time diversity available. This case is illustrated in Tab. Ill 
for L — 8, where it is shown that increasing the number of states results in a larger diversity. 
Note also that the optimum P-STC are able to achieve a diversity equal to that achieved by 
using the optimal generators for the AWGN channel, and that are not able to reach the diversity 
achievable on the RBFC with nL fading levels per codeword. This means that convolutional 
codes are more capable to collect time diversity than spatial diversity. 

Finally, we report in Tabs. IV, V, and VI the optimum generators obtained through the search 
algorithm for n = 2,3,4, respectively, with BPSK and QPSK modulations, and for different 
number of states. The corresponding performance factors are also reported. It is worth noting 
that, although the codes are not geometrically uniform, in most cases the code search based 
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on the ST-GTF with a fixed reference codeword leads to the same code as the search over all 
possible transmitted codewords, or to a code with similar performance. 

VIII. Conclusions 

In this paper we investigated the feasibility of a pragmatic approach to space-time codes, where 
common convolutional encoders and decoders are used, with suitably defined branch metrics. 

We extended to pragmatic space-time codes the concept of generalized transfer function for 
convolutional codes in block fading channels, that results in the possibility to rank different codes 
with an efficient algorithm, based on the asymptotic error probability union bound. A search 
methodology to obtain optimum generators for different fading rates has then been proposed. 

It has been shown that P-STC achieve better performance compared to previously known STC 
and that they are suitable for systems with different spectral efficiencies, number of antennas 
and fading rates, and are therefore a valuable choice both in terms of implementation complexity 
and performance. 
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Appendix 

An example of computation of the ST-GTF for P-STC is reported here. We consider a (2,1,2) 
2-PSTC used with m — 1 receiving antenna over a quasi-static BFC and obtained from a 1/2 
convolutional code with generators (1, 3) 8 and BPSK modulation format. These generators are the 
best for two-states CC in AWGN, with free distance 3; when used to build a P-STC we checked 
that this choice of generators produces the second best code in quasi-static fading channels, 
achieving diversity 2, i.e. the maximum available diversity. This code has been chosen since 
simple enough to allow the evaluation of time — ST-GTF, T M0 (D), using standard algorithm 
on the modified error state diagram. The two-states trellis is depicted in Fig. 7 (up-left). 

The associated possible output symbols [ci;c 2 ] and [g\]g2\ are in the set {X , Xi, X 2 , X 3 } 
with X = [-1; -1], Xi = [-1; 1], X 2 = [1; -1] and X 3 = [1; 1]. Thus, the matrix A(c,g) is in 
the set {a, b, c, d, e} with 

'°°) b=H c=H d=H eJ 4 - 4 

The error state diagram modified by splitting the all-zero state with different labeling 6 and 6 
is given in Fig. 7 (up-right). Thus, by following the steps in Sec.V-A we rewrite the error state 
diagram for quasi- static fading channel as in Fig. 7 (down-left) where 

/ D c D c \ ( D h D h \ ( D d D e 

<*(£>) = 1/2 /*(£>) = 1/2 h J 1(D) = 1/2 [ 

\D C D c j \D h D h j \D A D e 

The corresponding ST-GTF results in 

Q T S (D) = V T a(D) + Q T S (3(D) ^ j Q T S (D) = V T a(D) (I - ^D))' 1 
Ql(D) = Q T s(DHD) ^ 1 QJ(D) = U^a(D) (I - ^D))^ 7p) 



(33) 



giving 

i D c+d 
Tmo(D) = Ql(D)U = -[D* D c ] (I - P{D)Y l [D d D d ] T = — — , (34) 

which can also be expanded in 

T M0 (D) = D c+d (l + D b /{1 - D h )) = D c+d + £> c+d+b (l + D b /(1 - D h )) = ... 

= D c+d + D c+d+h + D c+d+2h + ... . (35) 
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Therefore, to obtain the ST-GTF of the code we need to compute the eigenvalues of matrices 
of the form c + d + ib with % > integer, and 



+00 



T[T M0 (D)] = J2di ll d^\ 



(36) 

6 ± while their 



i=0 

where Aji and \ i2 are the two eigenvalues of c + d + ib (e.g., A i = A 2 
product is simply \n\i 2 — 2i + 1. Then, we obtain 77 min = 2 and F miD (l) xs (JV/16)(1 + 1/3 + 
1/5 + . . .). It is interesting to note that both a(D) and (3(D) have equal elements, that e and 
e can be interchanged without altering the eigenvalues, and thus the ST-GTF does not depend 
on the particular reference sequence. Therefore this code is geometrically uniform even at the 
output of the channel. As check we can evaluate T (D) for the all-zero reference codeword 
c = Co by exploiting the error state diagram with scalar labels; to this aim it is sufficient to 
replace a(D) = D c , /3(D) = D h and 7(D) = D d , so obtaining 

D c+d 



T (D) =«(£>)(! -P(D))-^(D) 



(37) 



l-D b ' 

which is equal to T M0 (D) in (34). 

To evaluate the ST-GTF of the P-STC in a BFC with L = 2 and periodical interleaving, i.e., 
A = (DiD 2 D 1 D 2 . . .), we can extend the error state diagram by replacing the state S with two 
states Si and S 2 , thus obtaining transitions — > S 1 with output D\, Si — > S 2 with output D^, 
S 2 -> Si with output D\, Si -> 6 with output D% and S 2 -> 6 with output Df (see Fig. 7 
down-right). The ST-GTF is then 



Q Sl (Di, D 2 ) = Dt + Q S2 (Di, D 2 )D\ 
Q( h (D u D 2 ) = Q Sl (D u D 2 )D* 
Q 6 (Di, D 2 ) = Q Sl (Di, D 2 )D$ + Qs 2 (D u D 2 )Df 

and after expansion 



Q8 1 (D U D 2 ) = T= ^ B§ 
Qs 2 (Di,D 2 ) 



D c iD b 2 
ic nbi n c + d r>b 



[Q (Di,D 2 ) = ^±^ 



(38) 



To^, 7J 2 ) = T M0 (Di, D 2 ) = D\D d 2 + D c + d D* + D{ +h D d+h 



D c+d+h D 2h + ^ (39) 



in which only D\D A has two terms with rank 1. Therefore T[D\Df\ = d\id\i implying r] min = 2 
and F min (l) w JV/16. 
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Fig. 1. Equivalent low-pass scheme for space-time codes. 
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Fig. 2. Equivalent low-pass scheme for the proposed pragmatic space-time codes. 
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input bits 




Fig. 3. Example of 1 bps/ Hz pragmatic space-time encoder for n — 2 transmitting antennas and BPSK modulation, obtained 
with a rate 1/2 convolutional encoder with 4 states and generator polynomials (5, 7)s. 




Fig. 4. Trellis for the pragmatic space-time code of Fig. 3. On each branch the first is the input bit, followed by the two output 
antipodal symbols. 
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Fig. 5. Example of 2 bps/Hz pragmatic space-time encoder for n — 2 transmitting antennas and QPSK modulation, obtained 
with a rate 2/4 convolutional encoder with 4 states and generator polynomials (06, 13, 11, 16)s. 
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output bits 



Demod. 


' m 


j(-i) 






branch metric 



Er=i 1^ - VK(hi, a ci + h 2iS c 2 ) 



Fig. 6. Receiver structure for the proposed pragmatic space-time codes. In the figure the Viterbi decoder is the usual for 
the convolutional code adopted in transmission, with the only change that the metric on a generic branch is, for n = 2, 
J2T=i ~ [h^ a ci + ^2^C2 ) 2 , being £1,62 the two symbols associated to the branch. 
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Fig. 8. FER vs. SNR for P-STC obtained with the rate 1/2, 64 states convolutional encoder with generators (133, 171)s, 
1 bps /Hz BPSK, n = 2 transmitting antennas and m = 1, 2, 4 receiving antennas in BFC for various L. 
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SNR (dB) 

Fig. 9. FER vs. SNR for P-STC obtained with 1/4 convolutional encoder, 64 states, BPSK, 1 bps /Hz, 4 transmitting antennas 
and 1,2,4 receiving antennas in quasi-static Rayleigh fading channel. 
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TABLE I 

Comparison of rate 2/4 P-STC with QPSK, n = m = 2, n = 2 on BFC with L = 1, S 6 /iVo = 15<±B and N = 130. 

The performance factor is truncated with 8 h = 9. The first two codes are the best produced by the 
search (there are 12 first-class codes with almost the same behavior, and 48 second-class codes). the 

third code is the code proposed in [8]. 



Generators 


df 




F min (2)/iV 


Fmi„(c ,2)/iV 


p 

1 woo 


FER 


(06,13,11,16)8 
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0.076 


0.048 


3.5 10~ 4 


9.3 10~ 5 


(05,11,06,16) 8 
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4 


0.073 


0.092 


3.4 10~ 4 


1.0 10~ 4 


(01,02,04, 10) 8 
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0.125 


0.125 


5.7 10~ 4 


2.4 10~ 4 




Fig. 10. Exhaustive search for MIMO(2,2) P-STC in terms of FER: QPSK, quasi-static fading channel (L = 1), E b /N = 9dB. 
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01 23456789 10 

SNR (dB) 

Fig. 11. Comparison of our 4 states QPSK MIMO(2,2), 2bps/ Hz, P-STC (continuous line) and previously known space-time 
codes. TSC: [8], BBH: [12], and YB: [14], quasi-static fading channel. 
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Fig. 12. Comparison of our 4 states QPSK MIMO(2,2), 2bps/ Hz, P-STC and the STC in [8] (TSC) for different fading rates. 
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TABLE II 

Optimum rate 1/2 P-STC with BPSK, n = 2, m = 1 on BFC with L = 1. Parameters with superscript (1) refer 
to the codes obtained with best convolutional codes for AWGN channel. The performance factor is 
truncated with 8h = 2d^ — 1. Error performance refers to E b /N = 2QdB and N = 130. 
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TABLE III 

Optimum rate 1/2 P-STC with BPSK, n = 2, m = 1 on BFC with L = 8. Parameters with superscript (1) refer 
to the codes obtained with best convolutional codes for AWGN CHANNEL. The performance factor is 

TRUNCATED WITH Sh = 2e^ 1) — 1. fj^£ FC) IS THE DIVERSITY ACHIEVABLE ON THE RBFC WITH nL FADING BLOCKS 

[24]; THE VALUE OF THE SINGLETON BOUND IS 9. 
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TABLE IV 

Optimum P-STC for a system with n = 2, m = 1 on BFC with L = 1. Symbol * indicates that the search 

BASED ON F min (c , m) LEADS TO THE SAME CODE AS THE FULL SEARCH. 
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TABLE V 

Optimum P-STC for a system with n = 3, m = 1 on BFC with L = 1. Symbol * indicates that the search 

BASED ON F min (c , m) LEADS TO THE SAME CODE AS THE FULL SEARCH. N.E.=NOT EVALUATED. 
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TABLE VI 

Optimum P-STC for a system with n = 4, m = 1 on BFC with L = 1. Symbol * indicates that the search 

BASED ON F min (c , m) LEADS TO THE SAME CODE AS THE FULL SEARCH. N.E.=NOT EVALUATED. 
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